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In studies of island nucleation and growth, the distribution of capture zones, essentially proximity 
cells, can give more insight than island-size distributions. In contrast to the complicated expressions, 
ad hoc or derived from rate equations, usually used, we find the capture-zone distribution can be 
described by a simple expression generalizing the Wigner surmise from random matrix theory that 
accounts for the distribution of spacings in a host of fluctuation phenomena. Furthermore, its 
single adjustable parameter can be simply related to the critical nucleus of growth models and 
the substrate dimensionality. We compare with extensive published kinetic Monte Carlo data and 
limited experimental data. A phenomenological theory sheds light on the result. 

PACS numbers: 68.35.-p,81.15.Aa,05.40.-a,05.10.Gg 



In the active field of statistical mechanics applied to 
materials, an important unsettled problem in morpho- 
logical evolution during epitaxial thin film growth [l| is 
the characterization of the statistical properties of nucle- 
ating islands. In particular, for over a decade the univer- 
sal scaling shape of the island-size distribution (ISD) has 
been investigated numerically with kinetic Monte Carlo 
(kMC) simulations, but analytical evaluation has proved 
elusive. Only rate equations 0, 0, 0| or complicated (of- 
ten implicit) expressions [H, [y] have been proposed. The 
ISD is an important tool for experimentalists, since sim- 
ulations have shown it to be a unique function of the 
size i of the critical nucleus (see below) , a quantity that 
describes the largest unstable cluster. 

A decade ago Blackman and Mulheran 0, @] proposed 
subordinating the ISD to the distribution of areas of 
Voronoi polygons (proximity cells) built around the nu- 
cleation centers. Once an island is nucleated, it efficiently 
captures most of the adatoms diffusing within the cap- 
ture zone (CZ), a region roughly coinciding with the is- 
land's Voronoi polygon. This breakthrough, which in- 
vites analogies to distributions of quantum dots Q and 
to foams [9j, led to several part-numerical, part-analytic 
investigations 0, 0, 0, 0, d, Q that allowed prediction 
of the ISD for point islands with good accuracy, at the 
price of performing extensive kMC simulations or of solv- 
ing a system of several coupled, non-linear rate equations, 
which is computationally as taxing as kMC. For this rea- 
son, an empirical functional form, proposed in Ref. Q, 
which fits kMC results well, is still widely used to analyze 
data. 

In this Letter, we propose a different approach. We 
show that the generalized Wigner surmise (GWS) dis- 
tribution, a class of probability distribution functions 
rooted in random matrix theory (RMT) [l(3> [ll[ , yields 
an excellent quantitative description of the CZ size dis- 
tributions for all values of the critical nucleus size i in 
published simulations. Thus, this relatively mature sub- 



ject is vitalized and broadened by linkage to universal 
aspects of fluctuations. RMT experts will find remark- 
able that the signature exponent has atomistic meaning 
in these non-equilibrium systems. A phenomenological 
argument suggests the physical origins of the GWS here. 

RMT [l(J, EH has been successfully applied as a phe- 
nomenological description of statistical fluctuations in 
a large variety of physical systems, such as highly ex- 
cited energy levels of atomic nuclei, quantum chaos 
cross-correlations in financial data [13J, stepped crystal 
surfaces [3], and even times between buses in Cuernavaca 
and distances between parked cars [161 ] ! The last ex- 
ample is analogous to our study in that the RMT-derived 
formula accounts for the data notably better than the 
culmination of years of problem-specific theories. 

RMT applies only to matrices with special symme- 
tries, which constrains the applications to physical sys- 
tems that somehow reflect these symmetry properties. 
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FIG. 1: [Color online] Plots of the GWS of Eq. © P n (s), 
n=l,2,3,4, of relevance in this paper, indicated by long dashes 
alternating with n short dashes; also Ps/2(s), indicated by 
long dash, short dash, dot. The thin solid [blue] curve shows 
the Gamma distribution Il7(s), discussed later. 
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FIG. 2: a) Schematic for ID (vertical). Black rectangles cor- 
respond to ID islands. Horizontal lines mark the midpoints 
between the edges of two neighboring islands, with the capture 
zones (CZ) defined as the resulting proximity cells. An alter- 
native definition, implicit in the point-island approximation, 
uses the midpoint between the centers of islands, indicated 
by dashed lines and leading to Voronoi "cells". For islands 
[nearly] the same size, the two [nearly] coincide, b) 2D il- 
lustration of the islands (approximated as circular) and the 
Voronoi polygons that bound their CZ, from Ref. Q. 



The Wigner surmise Pp{s) 



Pp(s) = aps 13 exp(-bps 2 ) 



(1) 



(cf. Fig.[T]) provides a simple, excellent approximation for 
the distribution of spacings for such cases Hi 11 1. Here 
s is the fluctuating variable divided by its mean, f3 is the 
sole WS parameter [17[, and ap and bp are fixed by the 
normalization and the unit- mean conditions [l8| . 

Standard RMT 1^, 11 , 12} fixes attention on the values 
1, 2 and 4 of /3, corresponding to orthogonal, unitary, and 
symplectic matrices, respectively. The GWS posits that 
Eq. JT]) has physical relevance for general non-negative (3 
[19| . We show here that the CZ distribution is excellently 
described by the GWS with ( 2 /d)(i+l), where d=l,2 is the 
spatial dimension (see Fig. [2]) . The GWS also describes 
the distribution of terrace widths on stepped surfaces [3, 
[Hi} , where the step-repulsion strength determines (3. 

The explicit dependence of (3 on dimension is a novel 
feature of this study. Most other applications of RMT 
are either essentially one-dimensional or insensitive to 
d. While -Pi(s) describes well the nearest- neighbor spac- 
ings between randomly-distributed points on a plane [l2| , 
GWS fits of Voronoi tessellations of such points are not 
particularly good [2(| [2l[. For island nucleation, subtle 
correlations between nucleation centers lead to a distri- 
bution of tesselation cells (Fig. [2]) described by the GWS. 

Island nucleation is pictured as atoms deposited on a 
substrate (at rate F) and then diffusing on the surface 
at diffusion rate D (most properties depending only on 
D/F [![). When adatoms meet, they form bonds, whose 
lifetime depends on temperature T . At low enough T, 




FIG. 3: [Color online] CZ size distributions for critical nucleus 
size i — 1 and d = l. Symbols are from Fig. 12 of Ref. [e|, for 
various pairs of values of (D/F in units of 10 s , coverage in 
monolayers): □ (5, 0.11), (5, 1.19), A (5, 12.65), x (50, 
0.11), + (50, 1.19). The thin curve is the theory prediction of 
Ref. [g]. The thick [red] curve is the simpler P4(s). The inset 
shows similar results for distribution of gaps between point 
islands from Fig. 11 of Ref. [5J, with the added thick [red] 
curve giving Ps/2{s); its self-convolution is about Pi(s), the 
ISD, as discussed near the end. 



bonding is virtually irreversible, so that an adatom pair 
is a stable — and immobile island, which grows only by 
capturing other adatoms. A single adatom is then called 
a critical nucleus] equivalently, the critical nucleus size is 
i = 1 at low T . At higher T a single bond will be broken 
before other adatoms can be captured, and the critical 
nucleus will be a larger cluster, whose size will depend 
on the surface lattice symmetry, generally i = 2 or 3 on 
a (111) or (100) surface, respectively [2LI22I]. 

We first test our approach on data computed by Black- 
man and Mulheran [5| with kMC simulations of the nu- 
cleation of point islands along a one-dimensional (ID) 
substrate (cf. Fig. [2^,). Since i—1 there, we predict that 
the CZ size distribution is a GWS function with /3=4. 
Fig. [3] shows the results of their simulations, along with 
fits with the Wigner surmise. Clearly Pa(s) yields an ex- 
cellent fit to the numerical results, better than the thin 
solid line, the result [1] of a statistical numerical calcu- 
lation replacing the solution of a complicated integro- 
differential equation. Thus, our expression is both more 
accurate and simpler than their theoretical result. 

Two-dimensional (2D) deposition, diffusion, and ag- 
gregation models have been extensively treated by many 
authors. Mulheran and Blackman [3] report kMC sim- 
ulations of growth of fractal islands (i — 1) and circular 
islands (i = 1 to 3). For the circular islands we find very 
good agreement between the data and the GWS using 
(3 = ( 2 /2)(* + l) [21], w fih the trend for increasing i well 
reproduced. Even better agreement is found between the 
GWS P i+ i(s) and Mulheran and Robbie's @ more re- 
cent kMC simulations of nucleation and growth of circu- 
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FIG. 4: [Color online] (a) Symbols are numerical data from 
Fig. 2b of Ref. [g|, giving the CZ size distribution for nucle- 
ation of islands with i — in 2D. The thick (blue) curve is 
Pi(s). (b) Same as panel (a), but symbols for i = 1 from 
Fig. 2d of Ref. [(J, and the thick curve is ^(s). In both 
panels the thin curve is the theory of Ref. @|. 



lar islands for i = and 1, as shown in Fig. QJi and 2)d, 
again superior to their numerical-analytical theory [6(. 
Popescu et al. [24j also report extensive kMC simulation 
data of irreversible nucleation (i = 1) of point, compact, 
and fractal islands, but do not compute CZ size distri- 
butions. Their rate-equation approach was designed to 
describe island sizes and capture numbers, so should not, 
and does not (23[, describe the CZ distribution well. 

To understand why the CZ distribution is well de- 
scribed by Pp(s) with j3 = ( 2 /rf)(« + l), we offer a phe- 
nomenological model. We draw on our recent demonstra- 
tion [25| that the GWS appears in the context of RMT 
as the mean-field solution of Dyson's Brownian motion 
model ll|, |l2| , based on a Coulomb gas of logarithmically 
interacting particles (2(| in a ID quadratic potential well. 
We argue that the CZ size distribution can be extracted 
from a Langevin equation for a fluctuating CZ size in a 
confining potential well created by two competing effects: 
1) The effective confining potential well should increase 
for small-size CZ: nucleation of a new, small island causes 
a CZ of finite (and not greatly different from the mean) 
size to appear (cf. Fig. 19 of [27|), so that a large force 
must prevent fluctuations of the CZ size towards vanish- 
ing small values 2) The neighboring CZ's also prevent 
the one under scrutiny from growing, exerting a sort of 
external pressure, which may be assumed to come from 
a quadratic potential. A noise term represents atoms in 
a CZ attaching to other than the proximate island. 

To compute the repulsion, we analyze quantitatively 
the nucleation of new islands, following Ref. [28|. If N 



is the stable island density, n the adatom density, D the 
adatom diffusion constant, a the capture coefficient of an 
island, and Ni the density of critical nucl ei ( islands with 
i atoms), the nucleation rate N in 2D is [2a ] 

N = cmNi « Dn l+1 . (2) 

In the rightmost expression we have used a ~ D [28| and 
the Walton relation TV* w ri 1 [291] . In 2D, n satisfies 



F - anN « 



F/{<tN) w Ft 2 ID (3) 



where 1/N w I 2 = Dt is the squared diffusion length 
of adatoms before capture by an island in lifetime t. In 
ID, Eq. ([2]) still holds, but accounting for the properties 
of a random walk in ID makes the capture coefficient a 
dependent on I. Indeed, from jjl 30|] we have 



N = a lv> nNi 



Ni(D r ) 1/2 « jnNi 



D 



i+l 



(4) 



whence a%r> = D/l, Since n oc £ 2 , regardless of d, and 
1/N w £ in ID, Eqs. © and flU) can, for d=l,2, be written 



TV 



■.i+l 



D/l 



2-d 



(5) 



Taking the result for N/a in Eq. (JSJ) as a multiplicity, 
we find an effective entropy S = ks ln(n J+1 ). Identifying 
s = l d as the "area" in e?=l,2, and recalling n oc I 2 



d(Z/k B ) 
ds 



K 



i 2 ld)(i + l) 



(6) 



where if is a kinetic coefficient. Fluctuating repulsion 
(with strength B) from the neighboring CZs yields a sec- 
ond contribution (s)2 = —KBs + 77, where rj arises from 
the random component of the external pressure. We ar- 
rive at the Langevin equation 



s = K[(y d )(i + l)/s-Bs}+r,. 



(7) 



As we show in Ref. [25|, the stationary solution to the 
Fokker-Planck equation corresponding to Eq. (O is just 
the GWS P (s), with /?=( 2 Mi+l). 

Mulheran and Blackman [7|, |31[ proposed the Gamma 
distribution as an empirical description of general 
Voronoi tessellations, particularly of the CZ size distri- 
bution. With unit mean enforced, it has the form 32] 



n Q (x) = [a a /T(a)} x"- 1 exp(-ax). 



(8) 



No relation was established between the parameter a and 
any nucleation property, and in accounting for a range 
of i, a single value of a associated with random depo- 
sition is used 0]. The GWS and Gamma distributions 
are qualitatively similar, and, for 1 < (3 < 4, a is roughly 
2/3 + OiQ, where a is an offset of order one (cf. Ps(s) 
and n 7 (s) in Fig. {!])); the value of a depends on what 
property of the two distributions are equated [33| . How- 
ever, the slower decay of the II Q (s) leads to considerably 
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greater skewness, with a distinctly greater shift of the 
peak to smaller s. Like Pp(s), n Q (x) approaches a Gaus- 
sian for large a. Trying to distinguish the two forms in 
the large-s tail is problematic since the small values lead 
to large fractional uncertainty. Very recently H a (s) has 
been used as a tool for analyzing experimental CZ dis- 
tributions @, [H. Trial fits of the data with the GWS 
form are generally at least as good. Amar et al.'s pop- 
ular rate-equation-derived expression for ISD's, noted at 
the outset, is fi(s) cx s l exp(— idiS 1 ^), £>1, where aj is 
a complicated constant @. By construction, it peaks at 
s = l. While not designed for CZ distributions, /»(s) has 
been tried as an alternative to II Q (s) for quantum dots, 
with neither being fully satisfactory [8(. 

Several extensions and challenges present themselves: 
1) We hope experimentalists [1, 13] will refit their data 
for CZ distributions with the GWS. When Voronoi tessel- 
lation is used more generally, e.g. in studying biological 
systems [35j |. the resulting size histograms should be an- 
alyzed using the GWS and the alternatives. 2) Fig. [2^ 
shows that when there is a wide range of island sizes, us- 
ing a strict Voronoi construction rather than physically 
more appropriate edge cells [J, H?| leads to a narrower 
distribution (hence a higher deduced 0). E.g., in fitting 
numerical data for the irreversible point-island nucleation 
(i=l) in Fig. 3b of Ref. [53], we find much better agree- 
ment with Pi(s) than the expected ^(s). (Other sub- 
tleties complicate this case; compact islands occur only 
at very low coverage (<0.01) [27j.) 3) Our predictions 
should be tested in cases with large i. 4) For d >2, 
random walks are not recurrent, with the upshot that 
/3=i+l as for d=2. This can be tested with growth simu- 
lations in d=3 and d—A Q. 5)Ind=lit is straightforward 
to find the gap distribution between islands, as done in 
Ref.0. The CZ distribution is a self-convolution of the 
gap distribution 0]. In Fig. 1 we saw /3=4; the gap dis- 
tribution, then, should be well described by P 3 /2{s) (viz. 

P 4 (s) w 2 / 2s P 3/2 ( S ')P 3/2 (2 S ~ S ')ds'). The inset of Fig. 1 
corroborates this. Generalization to d=2 is unclear. 

In summary, as for spacings between parked cars [l6j |. 
the Wigner surmise provides a simple, universal expres- 
sion that accounts better for data than more complicated 
expressions developed over years of investigation. For our 
problem of the capture zone distribution in island nucle- 
ation, the exponent /3 of the generalized Wigner surmise 
P/3(s) provides information about the size i of the criti- 
cal nucleus and reflects the dimensionality d. Our phe- 
nomenological argument provides insight into the physi- 
cal origin of this behavior. Both features are significant 
advances beyond previous empirical analytic descriptions 
of the CZ size distribution (notably II Q (x)). The connec- 
tion to universal properties of fluctuations enhances the 
interest and importance of studies of CZ distributions 
and suggests many avenues for further investigations. 
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